Surface effects in doping a Mott insulator 
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The physics of doping a Mott insulator is investigated in the presence of a solid-vacuum interface. Using the 
embedding approach for dynamical mean field theory we show that the change in surface spectral evolution in 
a doped Mott insulator is driven by a combination of charge transfer effects and enhanced correlation effects. 
Approaching a Mott insulating phase from the metallic side, we show that a dead layer forms at the surface of 
the solid, where quasiparticle amplitudes are exponentially suppressed. Surface coiTelation and charge transfer 
effects can be strongly impacted by changes of the hopping integrals at the surface. 

PACS numbers: 71.38.-k, 71.30.+h, 73.20.-r, 71.38.Ht 



I. INTRODUCTION 

Recently, interest in the electronic properties of surfaces 
and interfaces of strongly correlated electron systems has 
grown due to the discovery of conducting interfaces, with 
a substantial high carrier mobility, between two insulating 
perovskites.' Theoretical studies of these systems^"^ con- 
cluded that the leakage of charge from one layer to the other 
explains the results. Thus, charge transfer between materials 
in heterostructures can be used to stabilize interface states that 
otherwise would be obtained only via the chemical doping of 
the parent oxide compound. On the other hand, comparing re- 
sults from measurements with different surface-sensitive pho- 
toemission, it has been shown that the spectral function can be 
strongly modified at the surface, compared with the bulk.^"^ 
For example, it was shown that the weight of the coherent 
part of the spectral function is reduced at the surface. ^'^ 

Interfaces obviously break translational invariance, and 
where present, the rotational symmetry of the bulk. Beyond 
this significant alteration in symmetry, there are important 
quantitative changes, such as a reduction of dispersion and 
changes of the crystal fields. Stress or strain is induced by 
the interface, which alters the distances and bonds between 
the ions. The lattice may even be structurally or electronically 
reconstructed. Due to these and other phenomena, the inter- 
face shifts and distorts the electronic states and energy levels 
at the ions, and modifies the entire electronic band. Since 
the electronic states of the individual ions and the screening 
are altered at interfaces, their correlations are strongly modi- 
fied as well.** These phenomena raise an important question: 
how does the surface/interface electronic phase differ from the 
bulk? 

Theoretical research on the effects of a surface or interface 
on strongly correlated electrons in inhomogeneous systems 
was initiated in large part by Potthoff and Nolting,^ who stud- 
ied the Hubbard model using an extension of dynamical mean 
field theory (DMFT)"^ appropriate to a layered system with 
a surface. Studies of solid-vacuum interfaces have unveiled 
the possibility of surface ferromagnetism, ' ' and explored the 
penetration depth of a bulk metallic phase into an otherwise 
insulating surface.'^''-' Borghi et al. have indeed shown the 
existence of a dead layer, due to an exponential penetration of 
metallic excitations.''* Yunoki et al}^ have studied the possi- 
bility of magnetic ordering and superconductivity at the in- 



terface between undoped high-Tc cuprates and manganites. 
Ishida et al. have studied the effect of interplanar Coulomb 
correlations on the electronic structure of strongly correlated 
heterostructures by applying the cluster DMFT method.'* 
Studies of solid-vacuum interfaces of the Holstein-like mod- 
els in an inhomogeneous system have shown that the polaron 
crossover takes place first on the surface. '^'''^ 

Most of this previous work has been done for the half- 
filled Hubbard model with particle-hole symmetry. However, 
at interfaces there is not only a change of the energy spec- 
trum and of the correlation, but due to charge transfer, the 
state occupancy also varies. In the present work, we focus on 
charge transfer phenomena in the framework of a Hubbard- 
like model Hamitonian. Especially, we are interested in a 
doping-driven Mott transition in the presence of a surface. 

The outline of this paper is as follows. In section II we in- 
troduce the model Hamiltonian, which is a semi-infinite Hub- 
bard model with layer dependent parameters. The correspond- 
ing results are presented and analyzed in section III. Finally, 
section IV is devoted to concluding remarks. 



II. THE MODEL AND METHOD 

The simplest model Hamiltonian exhibiting Mott physics'^ 
is the Hubbard model, in which the local repulsion corre- 
lates the electronic motion. At half-filling, this model cap- 
tures many properties of the pressure-driven Mott transition. 
The Hubbard model away from half-filling can be used to in- 
vestigate the density-driven Mott transition. Here we try to 
model some part of the phenomena occurring at the interfaces 
of the strongly correlated electron systems using the Hubbard 
model with non-uniform model parameters. The system is 
a three-dimensional, bipartite simple-cubic (sc) lattice with 
nearest-neighbor hopping only. The lattice is cut along a plane 
perpendicular to one of the coordinate axes, e.g. the z-axis 
[sc(OOl) surface]. The system is considered to be built up by 
two-dimensional layers parallel to the surface. Accordingly, 
the position vector to a particular site in the semi-infinite lat- 
tice is written as RsUe = i'i + Ra- Here Ra stands for the 
coordinate origin in the layer a, and the layer index runs from 
a = 1 for the topmost surface layer to infinity. The variable r j 
is the position vector with respect to a layer-dependent origin, 
and runs over the sites within the layer Each lattice site is then 
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labeled by indices i and a. In this notation, the Hamiltonian 
reads: 

H=- ^ Ua^jpc\^^Cj[3„ +'^Uan^a^nial, (1) 

{ia,jP)(j ia 

where Ciaa (^c|^cr) '■^^ destruction (creation) operator for 
electrons with spin a on site i of the a layer The electron den- 
sity on site ia is denoted riia, and tiajp is the hopping matrix 
element between two nearest-neighbor sites. Ua is the layer- 
dependent electron-electron on-site interaction. We use a uni- 
form U and fix the energy scale by setting i(iQj,g) = t = 1 
for a, /3 1. The intralayer hopping ta ji = ^ii or the hop- 
ping between the surface and the subsurface layer j2 = ti2 
can be modified (in, ti2 7^ t) to simulate different correlation 
strengths and/or relaxation processes at the surface. 

Our calculations are based on the embedding approach'^ '^ 
for dynamical mean field theory for a simple cubic lattice. 
This method is based on partitioning of the layered structure 
into a surface region which includes the first A'^ layers and 
an adjacent semi-infinite bulk region (substrate) to which it is 
coupled. Then the effect of the substrate on the surface region 
is described by an energy-dependent embedding potential. In 
constructing the embedding potential the self-energy of a bulk 
crystal corresponding to the substrate is required; this is ob- 
tained through a standard DMFT calculation. After construct- 
ing the embedding potential the self energy matrix of the sur- 
face region is determined using a layered DMFT calculation. 
In the framework of single-site DMFT, the self-energy matrix 
of the surface region is local [i.e., Y.ap{ioJn) ~ '^a{'i(^n)Sai3] 
and independent of wave vector. Using the converged self- 
energy matrix we can find quasiparticle weights and local oc- 
cupancies. 

In this study, the number of surface layers is chosen to be 
= 5 and we tested that this number provides converged 
results. Our impurity solver is exact diagonalization,^' where 
the bath is represented in terms of a finite number of levels, iig. 
The typical value we considered for the bath level is = 9 
and we tested that this number provides essentially converged 
results. All errors associated with the use of these numbers 
are smaller than the points on the graphs. 

III. RESULTS 

For a bulk lattice, in a phase with translational invariance, 
local occupations are site independent, (ua^) + (nai) = 
{ua) = n. On the contrary, for a semi-infinite lattice, the 
different local environment of the surface sites causes the lo- 
cal occupations near the surface to differ from the bulk filling. 
Charge transfer is present in the free-electron system-^ and 
we discuss the facts already known for this case first. The 
reduced coordination number and probably different hopping 
elements connecting at least one site in the surface layer of 
the semi-infinite system leads to a narrowing of the surface 
local density of states (LDOS), pfli{E), relative to the bulk, 

Pbuik i-^) (spectral weight is reduced at low and high energies 
but enhanced at intermediate energies). However, the surface 



LDOS and the bulk LDOS have the same band edges and their 
overall width is given by the free bulk DOS.-^ The band nar- 
rowing can be understood by referring to the moments of the 
LDOS, 

/OO 
E"'p^J{E)dE. (2) 
-00 

For a site in the bulk, the width of the DOS is given by 
^Pbaik = - (M(1'°))2 = 6*2, while for a site on the 

surface the result is ^p^^i = 4tf ^ + tl2-^* Clearly, if these 
hopping matrix elements are equal to those in the bulk, then 
this width is narrower at the surface than at the bulk. From this 
fact, it follows that in order to ensure thermodynamic equilib- 
rium, i.e., a common chemical potential /x, one has na=i < n 
below and na=i > n above half-filling, since the surface elec- 
tron density 2 pf}^-^{E)dE is lesser (greater) for below 

(above) half filUng than the bulk density, 2 p^°^i^{E)dE. 
Thus the surface layer will have charge depletion (accumula- 
tion) below (above) half filling. 

These effects are all due to a reduced coordination num- 
ber and probably different hopping elements. In addition a 
strong electron-electron interaction can cause a quite compli- 
cated charge redistribution at the surface. In the following, we 
present results on charge transfer and its effects on spectral 
weight for a semi-infinite Hubbard model with uniform- and 
non-uniform hopping elements. 



A. Uniform parameters 

We start our discussion with the case of uniform parame- 
ters. Fig. 1 shows the charge transfer Ana=i = na=i — n 
between the surface layer (a ~ 1) and the bulk as a func- 
tion of the filling n at T = 0. The three curves are for (a) 
U ~ as just discussed, (b) U/W = 1, a value chosen to be 
below the critical Mott transition value Uc/W w L33,'^ and 
(c) U/W ~ 1.5, a value chosen to be above the critical Mott 
transition value Uc- 

As argued above, the case without interactions displays sur- 
face charge depletion for all electron densities below half fill- 
ing (and exactly the opposite for above half filling). For cases 
where there are interactions between the electrons, the situ- 
ation is clearest in the strong coupling limit, where, at least 
at half filling, there are two well separated Hubbard bands. 
Insofar as doping into the 71 < 1 regime involves only the 
lower Hubbard band (definitely an oversimplification) then 
we should see an electron density dependence over this en- 
tire electron density regime (0 < n < 1) similar to what is 
observed in the case without interactions, but over the entire 
density range (0 < n < 2). This is borne out by the results 
in Fig. 1 for non-zero interaction strength: a surface electron 
density depletion is observed for small to intermediate densi- 
ties, and an accumulation for densities from beyond quarter 
filling to near half filling is seen. Apparently for an interact- 
ing system the key effect of a narrowing of the surface local 
density of states holds quite generally, even in the presence 
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FIG. 1: (Color online) The charge transfer AnQ=i = nc,=i — n 
between the surface layer (a = 1) and the bulk as a function of the 
filling n for uniform t. 



of interactions, and even within the lower and upper Hubbard 
bands. 

Upon increasing the bulk charge density, the chemical po- 
tential moves from the lower half of the lower quasiparticle 
sub-band to the upper half of that same sub-band. In addition, 
the differences between the 7i„ and the bulk n diminish with 
increasing distance from the surface and for the third layer, the 
occupation is almost indistinguishable from the bulk occupa- 
tion (not shown). Finally, note that the charge transfer in the 
noninteracting limit is larger than in the interacting cases, so 
the expectation that increased interactions tend to reduce over- 
all the magnitude of charge transfer is supported by these re- 
sults. Indeed, a finite charge transfer implies an (unfavorable) 
increase of the average double occupancy compared with the 
homogeneous charge distribution; thus, a non-zero Hubbard 
interaction between conduction electrons will tend to suppress 
charge transfer However, due to separation of low- and high- 
energy states, charge transfer is not completely suppressed in 
the strong coupling limit, U ^ oo, even when the long range 
Coulomb interaction is present. 

Also as expected, at half-filling there is no charge transfer at 
all, because particle-hole symmetry on a bipartite lattice with 
nearest-neighbor hopping requires a homogeneous charge dis- 
tribution among the layers. 

Theoretical studies of the Hubbard model for bulk sys- 
tems, mainly using DMFT,'° have shown that in a bulk sys- 
tem at half-filling, the system is insulating for a sufficiently 
large Coulomb repulsion, U > Uc- For any other filling, 
the ground state is metallic. In the metal the charge com- 
pressibility K ~ {\/7T?){dn/dn) remains finite even when the 
insulator is approached from above or below half-filling, but 
the chemical potential changes discontinuously from yu = /i^ 
to /.J = /i+ as the filling goes from n = 1^ to ?i = 1+. 
For 1^1^ < n < fi^, the system is insulating and half-filled. 
Doping away from half-filling induces states inside the Mott- 
Hubbard gap}^-^^ Therefore the jump in chemical potential 
Afj, = fi^ — for infinitesimal doping is strictly less than 
the single-particle gap of the insulator 

Another feature of the DMFT treatment of the Mott tran- 
sition is that both the interaction-driven transition at half- 
filling'" and the density-driven transition^*^"^^ are character- 
ized by the coexistence of two solutions in some parameter 
regimes. In the latter case, there is a region of chemical po- 



tentials, delimited by two curves fJ.~i{U) and n~2{U) for hole 
doping, in which a metallic solution with n ^ 1 and an in- 
sulating one with n = 1 coexist. Naturally, the stable so- 
lution is the one that minimizes the grand-canonical free en- 
ergy. At T = the metal is stable in the whole coexistence 
region and is continuously connected to the insulator At fi- 
nite temperature the energetic balance is more complicated, 
and a first-order transition occurs along an intermediate line 
between jici and fic2, leading to phase separation between the 
two phases (formation of a phase mixture with a fraction Xmet 
of the metal and 1 — Xmet of the insulator). 

Is it possible to drive a Mott transition on the surface at a 
critical chemical potential different from that required for the 
bulk? To address this question, we focus on the behaviour of 
the electron density n as a function of the chemical potential 
/i and on the density-driven Mott transition. A plot of the 
occupation as a function of the chemical potential, shown in 
the left panels of Fig. 2 for two different values of U, clearly 
shows the Mott transition into an incompressible state, as a 
plateau near /i = U/2, when the Hubbard U is greater than 
the critical coupling Uc. When U is lower than the critical 
value (upper left panel) then, while the filling still approaches 
unity at /i = [//2, there is clearly no plateau present. 

As can be seen in the lower left panel, for chemical potential 
near the critical value (where becomes unity), the surface 
layer approaches an insulating state slightly faster than the 
bulk. Nonetheless, the electron density eventually approaches 
unity at the same chemical potential that governs the bulk den- 
sity. The net result is that the surface tracks the bulk as far as 
the metal-insulator transition is concerned, and a metal in the 
bulk implies a metal at the surface. 

To examine this transition in more detail, the right pan- 
els of Fig. 2 show the quasiparticle weight for the a*'' layer, 
Zq = (1 — dT,a{i^)/duj\cu=Q) ^ CSaii^) IS the self-energy for 
layer a) as a function of occupation for U/W = 1.0 and 
U/W = 1.5. This figure shows the existence of a so called 
dead layer at the surface characterized by an exponentially de- 
caying quasi-particle weight near half-filling (lower panel). In 
the upper panel, the quasiparticle weight in the bulk does not 
reach zero, indicating that the U is below the critical value. 
Smaller z at the surface at small doping means the coherent 
peak near the Fermi energy is narrower at the surface and 
the incoherent Hubbard band is more pronounced than in the 
bulk, in agreement with experiment. Such a trend has been 
reported in studies of the correlation-driven metal-insulator 
phase transition at half-filling. ' ' However, the reduction of the 
surface quasiparticle weight seen here is much weaker than in 
the experiments. 

At low electron density the quasiparticle weight for the sur- 
face layer is greater than for the bulk (not shown). In this re- 
gion the depletion of charge from the surface to the bulk (see 
Fig. 1) causes a reduction of correlation effects, which leads 
to a more robust quasiparticle at the surface. 

To illustrate to what extent charge transfer is responsible 
for the differences at the surface vs. the bulk, we try to elim- 
inate the effect of enhanced correlation at the surface. In the 
left panel of Fig. 3 we show the bulk and surface quasi parti- 
cle weights for U/W — 1.5 as a function of electron density 
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FIG. 2: (Color online) Left panels: Occupation n vs (/i — f/ /2) jW for two different values of interaction strength U jW . The data shown are 
for the bulk system and the surface layer with uniform t. The plateau in the occupation at n = 1 signals the onset of the incompressible Mott 
state. Right panels: Surface layer quasiparticle weight z\ for uniform t and the bulk quasiparticle weight as a function of occupation. The 
solid line shows z for the bulk calculation. 



very close to half-filling. This figure also include several bulk 
calculations for larger values of U . If the enhanced correla- 
tion effect is responsible for the differences at the surface, we 
should recover the original value of the surface z in a bulk cal- 
culation with a larger value of U . As is seen from the figure, 
even for a very large value of [/, the bulk calculation is un- 
able to reproduce the original surface result. Therefore charge 
transfer effects play an important role in a doped Mott insu- 
lator As another test for this statement, in the right panel of 
Fig. 3 we use U jW = 2.0 in the bulk, and then show several 
semi-infinite calculations with reduced surface correlation IJ \ . 
This figure also shows that upon reducing the value of C/i it 
is not possible to have the same quasiparticle weights at the 
surface as in the bulk calculation. Thus, the charge transfer 
effect is more important than enhanced correlation effects in 
a doped Mott insulator This is expected because the interac- 
tion strength is already large enough to give a Mott insulator 
at half-filling, so enhancing that value would not lead to any 
significant change. Indeed, the occupation distance from half- 
filling is more important here. This quantity tunes with the 
charge transfer effect. 

B. Modified surface parameters 

In the previous section, we mentioned that the reduction of 
the quasiparticle weight at the surface is smaller when com- 
pared to the reduction found in experiments. Here we show 



that a modified hopping at the surface causes a larger charge 
transfer and then the reduction in surface quasiparticle weight 
is greatly amplified. In Fig. 4 we show the charge transfer 
for the case where the hopping from a site on the surface to 
one on the layer beneath is lower, = 0.5, and the case 
where the intralayer hopping on the surface layer is reduced 
with respect to that in the bulk, t\\ =0.5. These figures show 
that reduced hopping matrix elements on the surface sites en- 
hance the amount of charge transfer In the case of t\\ = 0.5, 
the value of charge transfer is sufficiently large that including 
the long-range Coulomb interaction will have some effect in 
reducing the total amount. For simplicity, however, we ex- 
clude this point here; we believe that including the long-range 
Coulomb interaction will lessen some of the effects described 
below, but the basic results will still hold. 

Fig. 5 shows the electron density as a function of the chem- 
ical potential (left panels) and quasiparticle weights as a func- 
tion of electron density (right panels) for reduced hopping pa- 
rameters near the surface. As can be seen in Fig. 5 the ten- 
dency toward an insulating state is certainly enhanced as we 
lower the hopping parameters near the surface. For reduced 
parameters, particularly for reduced tn, the plateau region is 
considerably extended (left panels). Close examination, how- 
ever, verifies that even in the case of non-uniform parameters, 
there is only one transition for the entire sample. Also, results 
for quasiparticle weights (right panels) show that the reduc- 
tion in z is considerably more pronounced in the case of re- 
duced parameters and is almost negligible for a wider range of 
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FIG. 4: (Color online) The charge transfer Ana=i — 7ia=i — n 
between the surface layer (a — 1) and the bulk as a function of the 
filling n for ti2/t = 0.5 and tn/t = 0.5 at T = 0.0. 



FIG. 3: (Color online) Plot of the quasiparticle weight, z vs. elec- 
tron density, n, near half filling for U/W = 1.5 (top panel), and 
U/W = 2.0 (bottom panel). The bulk result is given by the solid 
(black) curve. The corresponding result for the surface layer is given 
by the red (dashed) curve with open square symbols; the value of z 
at the surface is significantly reduced from that in the bulk. Since the 
coordination number is lower on the surface then the surface elec- 
trons effectively experience a large Coulomb repulsion. If this is the 
primary reason for the reduced value of z, then a bulk calculation 
with a larger value of U should mimic the original result at the sur- 
face. This is clearly not the case, even when we double the value of 
U. In the right panel we instead lower the value of Ui (on the sur- 
face) in an attempt to restore a value of z equal to that in the bulk; 
this is also unsuccessful, indicating that charge transfer effects are 
an important ingredient for the reduced quasiparticle weight on the 
surface. 



doping. It is a matter for experiment to decide if such changes 
in surface parameters can explain the data in a consistent man- 
ner. 

To further illustrate to what extent charge transfer is respon- 
sible for the differences at the surface vs. the bulk, even in 
cases where the surface parameters can vary from those in the 
bulk, we reshow in Fig. 6 the results for z for U/W = 1.5 
and tu/t = 0.5 for both the bulk (open squares joined by 
solid (black) lines) and the surface (filled red squares joined 
by dashed (red) lines). There are clear and significant differ- 
ences, particularly at low {n w 0) and high (n « 1) electron 
concentration; however, to account for charge transfer effects 
we use the information in Fig. 4 and show the surface quasi- 
particle weight on the same plot as a function of ni; this is 
shown by open circles joined by (green) short-dashed lines. 
The trend as a function of electron concentration now is qual- 
itatively similar to that of the bulk. However, one can improve 
upon this even further, since, in reality the electrons on the sur- 



face experience a stronger relative Coulomb repulsion due to 
the reduced intralayer hopping parameter. For this reason we 
include a calculation in the bulk for U/W ~ 3.0 (since the 
intralayer hopping in this layer is half what it is elsewhere); 
this is shown with open triangles joined by (pink) dot-dashed 
lines, and these results are in excellent quantitative agreement 
with the renormalized surface quasiparticle weights. Thus, 
we find that essentially the entire change of the quasiparticle 
weight at the surface can be accounted for once the charge 
transfer effects are taken into consideration. 



IV. CONCLUDING REMARKS 

In conclusion we have investigated the surface electronic 
properties of a doped Mott insulator by applying the embed- 
ding approach for DMFT. As a model system, we considered a 
semi-infinite Hubbard model. We have shown that the reduced 
coordination number and probably reduced hopping elements 
at the surface causes an accumulation of charge at the surface 
in the strong coupling limit, that makes the surface approach 
the condition of half-filling more quickly than the bulk. The 
accumulation of charge at the surface, in turn, makes the cor- 
relation effects more effective. In particular, we have demon- 
strated that the reduction of the quasiparticle weight detected 
by surface sensitive photoemission experiments of a doped 
Mott insulator are caused by both charge transfer and en- 
hanced correlation effects at the surface. The expected modifi- 
cation of the intra-layer hopping at the surface and inter-layer 
hopping between the surface and the subsurface layer lead to 
more significant differences between the surface and the bulk. 
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FIG. 5: (Color online) Left panels: Occupation n vs (/i — [/ /2)/W for two different values of interaction strength U /W . The data shown 
are for the bulk system and the surface layer with t\2/t = 0.5 and tw/t = 0.5. The plateau in the occupation at n = 1 signals the onset of 
the incompressible Mott state. Right panels: Surface layer quasiparticle weight zi for non-uniform t and the bulk quasiparticle weight as a 
function of occupation. The solid line shows z for the bulk calculation. 
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FIG. 6: (Color online) The quasiparticle weight is plotted vs electron 
density. The open squares joined by solid (black) lines denote the 
bulk weight as a function of the bulk electron density; this shows a 
characteristic decrease from n ~ to half-filling, where an insulat- 
ing state is present. The surface quasiparticle weight, z\ (filled red 
squares joined by dashed (red) lines) is shown vs. the same bulk elec- 
tron density, and shows considerable variation from the bulk. The 
same quantity, but plotted vs. the surface electron density, is shown 
by open circles joined by (green) short-dashed lines. This result is 
in qualitative agreement with the bulk result. An improvement is 
achieved by accounting for the effective higher Coulomb repulsion, 
as described in the text; this is shown with open triangles joined by 
(pink) dot-dashed lines; this result is in excellent quantitative agree- 
ment with the renormalized surface quasiparticle weights. We have 
used t\i/t — 0.5 for the calculations shown in this plot. 
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